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We present a measurement of the relative branching fraction, Rj i^, of B® —¥ J/^>/o(980), with 
/o(980) -> tt+tt", to the process B° J/ipij), with <f> -> K+K~. The J/^/o (980) final state 
corresponds to a CP-odd eigenstate of Bg that could be of interest in future studies of CP violation. 
Using 8 fb _1 of data recorded with the DO detector at the Fermilab Tevatron Collider, we find Rf /4, 
= 0.275 ± 0.041 (stat) ± 0.061 (syst). 

PACS numbers: 13.25.Hw, 14.40.Nd 



The CP-violatingphase angle in B® mixing, (j>i^^, 
has been measured [lj-|3| using B® — > J/yj(j) decays, and 
is statistically consistent with that predicted by the stan- 
dard model (SM) [|L Ignoring possible ambiguities in its 
hadronic structure 5[, the weak phase angle <^ css mea- 
sured in B s — > J/-0/o(98O) decay should be equal to the 
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Jli>4> 
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SM 
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where 8? M is the SM 



angle in the unitarity triangle for the B° s system that is 
analogous to the well-known angle (3 in the B® system, 
and 4>^ p is any additional phase arising from new physics 
in B® mixing. Measuring this phase angle through vari- 
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ous decay modes could help reduce its uncertainty. In 
particular, the decay products of B° s — > J/ifxj) are in 
an indefinite CP state, requiring CP-even and CP-odd 
components to be extracted through a time-dependent 
angular analysis. In contrast, the decay products in 
— > J/-0/o(98O) are in a CP-odd eigenstate, which can 
provide a more direct measurement of <f> c s css relative to 
B® —> J r /ip4> [6]. The precision of a 4> c s css measurement in 
the /o(980) channel is expected to be poorer than in the 4> 
channel because of the smaller decay branching ratio for 
this process. However, a complementary method of anal- 
ysis provides different systematic uncertainties, as well as 
an important cross check of the result from B® — > J / 

In this Article, we present a measurement of the rela- 
tive branching fraction (-R/ m) which is, based on theo- 
retical estimates, expected to be significant (7- 11|: 



R 



fo/'P 



_ r{B° 8 ^ J/^/ (980),/q(980) 



T(B° -> J/i/j(j), 4> -> K+K-) 



0.20- 0.40. 
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The LHCb Collaboration has reported [12j a first 

measurement of iZ/ M = 0.252lg Q32 ( stat )^o 033 ( s y s t)- 
The Belle Collaboration has measured the branching 
fraction B(B° -> J/ipf (980), /o(980) ir+ir-) = 
[1.16^;?i(stat)t o: i 5 7 (syst)l°;?6(7V B «„ i3 (, ) )] x 10" 4 0, 

where iV («)=(„) is the number of i?s ; _Bg ; pairs in the 
sample. The CDF Collaboration has also measured the 
relative branching fraction and finds Rf u = 0.257 ± 
0.020 (stat) ± 0.014 (syst) 0. We report a new mea- 
surement of the relative branching fraction using data 
collected with the DO detector at the Fermilab Tevatron 
Collider. 

To determine an absolute branching fraction requires 
an excellent understanding of efficiencies, other related 
branching fractions, cross sections, and integrated lumi- 
nosity. However, by measuring a relative branching frac- 
tion, terms common to both the B° — > J /tpf (980) and 
the Bg — > J/^4> decays cancel, giving: 
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N B o_ 



So/<t> 



■J I -0/o (980) 



trcco 

J /-4>f (980) 

<-rcco 
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Hence, just the yields N B o_ 



J/V/o(980) 



and N B o_ 



and 

their detection efficiencies, efcco J ^^ and ef c B c^ J ^ fo< " 980 \ 
are needed to measure Rf /<p- 

The DO detector is described in Ref. [TBI, and only 
those components that directly affect this measurement 
are discussed below. The tracking system consists of 
a silicon microstrip tracker (SMT) and a central fiber 
tracker (CFT), both located within a 1.9 T supercon- 
ducting solenoid magnet. The SMT has approximately 
800,000 individual strips, with typical pitch of 50 — 80 
/im, and a design optimized for tracking and vertex- 
ing capability within the pseudorapidity range \r/\ < 3 
|16| . The system has a six-barrel longitudinal structure, 
each barrel having four layers arranged axially around 
the beam pipe, interspersed with 16 radial disks. The 
CFT has eight thin coaxial barrels, each supporting two 
doublets of overlapping scintillating fibers of 0.835 mm 
diameter. One doublet is parallel to the collision axis, 
and the others alternate by ±3° relative to that axis. 
The muon system resides beyond a calorimeter that sur- 
rounds the inner tracking detectors, and consists of one 
layer of tracking detectors and scintillation trigger coun- 
ters before 1.8 T toroids, followed by two similar layers 
after the toroids. 

Approximately 8 fb _1 of integrated luminosity is used 
in this analysis. The data are divided into four time pe- 
riods, corresponding to different detector configurations 
and instantaneous luminosities, called Run Ha (1.4 fb _1 ), 
Run Ilbl (1.4 ftr 1 ), Run IIb2 (3.3 fb" 1 ), and Run IIb3 
(2.1 fb- 1 ). 

We search for B® — > J/i/)/o(980) candidates using the 
decay mode J/ip — > Events are collected using 

a mixture of single and dimuon triggers which have a 



similar trigger efficiency for both B° s — > J/-0/o(98O) and 
B° s —> J/tp(f>. Muon candidates must have transverse 
momentum pt > 1.5 GeV and be detected in the muon 
chambers within the toroidal magnet. In addition, each 
muon track must be associated with a track reconstructed 
by the CFT, and have at least one SMT hit. The J/tp 
candidates are formed from two oppositely charged muon 
candidates emanating from a common vertex, and have 
at least one of the muon candidates detected outside the 
toroidal magnet. 

All reconstructed tracks not associated with muons 
forming a J ftp candidate are considered in the recon- 
struction of /o(980) and <j) candidates. Since the DO 
detector has limited ability to separate kaons from pi- 
ons, tracks are assigned the pion mass when searching 
for B° s — > J/?/>/o(980) and the kaon mass when search- 
ing for B® — > J/ip<p, Charged tracks are required to 
have at least two CFT hits, at least two SMT hits, a 
total of at least eight SMT and CFT hits, and a mini- 
mum pt of 300 MeV. Any two oppositely charged tracks 
that have one track with transverse momentum px > 1-4 
GeV, an invariant mass 0.7 GeV < M^+ v - < 1.2 GeV 
or 1.0 GeV < M K + K - < 1.05 GeV, and are consistent 
with originating from a common vertex, are considered as 
/o(980) and 4> candidates, respectively. The /u + /i~7r+7r~ 
(/j + (x~~ K + K~) candidates are required to form a com- 
mon vertex and have an invariant mass between 5.0 and 
5.8 GeV. The invariant mass requirements on M 7T + 7r - and 
M K + K - prevent the two tracks to be considered as can- 
didates for both /o(980) and (f> interpretations. 

The final data sample is formed by applying the fol- 
lowing additional requirements to further reduce back- 
grounds. The /o(980) and cj) candidates must have pr 
> 1.6 GeV with 0.91 GeV < M OT+7r - < 1.05 GeV and 

1.01 GeV < M K+K - < 1.03 GeV. The B° s candidates 
are required to have pr > 5 GeV, 2.9 GeV < M M + M - < 

3.2 GeV, and have a proper decay length with a signifi- 
cance of greater than 5 standard deviations (sd). 

The proper decay length, defined as L xy ■ (M b o/pt), 
where pt is the transverse momentum of the B® , M R o is 
the PDG value of the mass of the B® [17], and L xy [3] is 
the transverse distance between the primary pp interac- 
tion vertex and the four-track vertex of the B® candidate, 
is calculated for candidate primary vertices that use the 
transverse beamspot as a constraint. If there is more than 
one such vertex in an event, the primary vertex nearest 
in the transverse plane to the J/ip candidate is chosen 
for this analysis. 

A final selection is based on two Boosted Decision Tree 
[H, [H (BDT) discriminants. We use the Monte Carlo 
(MC) pythia progra m 12 ll to generate B® events and 
the EVTGEN program [22[ to simulate their decay. MC 
signal and background samples are used to train a BDT 
and to form discriminant output values for each event. 
The expected background is primarily due to two sources: 
prompt background that is defined as directly produced 
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J/ip mesons accompanied by tracks from hadronization, 
and non-prompt, or inclusive B — > J/tp + X decays where 
the J/ip meson arises from a 6-hadron decay accompanied 
by tracks from hadronization. Two MC background sam- 
ples are therefore generated with pythia: a sample of di- 
rectly produced J / ip prompt events and an inclusive sam- 
ple of B° s for all decay processes B° J/tp + X. A MC 
signal sample of J/-0/o(98O) events is used to train both 
BDTs. Thirty input variables are used in the BDT, in- 
cluding the momenta of final-state objects, vertex-quality 
requirements, B° s isolation, and decay angles. Six BDT 
isolation variables are used in this BDT, representing 
different choices for the size and which tracks are in- 
cluded in the cone of isolation. The BDT selections for 
both prompt and inclusive training are defined with a re- 
quirement on the BDT output value which provides large 
S/\/B, while keeping the signal yields high, where S and 
B are the number of signal and background events. 

The invariant masses of /o(980) and B s — > J/ipfo(980) 
candidates, following BDT selections are shown in Figs.Q] 
and[2J respectively. An unbinned likelihood fit is used to 
determine the yield of signal in each sample. The /o(980) 
has a large width (T3] and a mass just below the KK 
threshold. This affects the line shape, which is not a 
simple Breit-Wigner form, particularly at large mass val- 
ues. The 7r + 7r~ mass distribution is therefore fitted using 
the functional form of Ref. [23j |. which takes account of 
the opening of the KK threshold, and is convoluted with 
a Gaussian resolution function with a sd of 15 MeV. The 
line shape determined by fitting the /o(980) in MC, us- 
ing a second-degree polynomial for the background is also 
used to fit the data. Candidates for B® — > J/-0/o(98O) 
are defined by the ir + ir~ invariant mass window 0.91 < 
M v+W - < 1.05 GeV. The B° s -> J/<0/ o (98O) mass distri- 
bution is fitted to a Gaussian signal, with a background 
function consisting of a second-degree polynomial and a 
Gaussian at lower invariant mass to take into account 
partially reconstructed B decays. The unbinned likeli- 
hood fit is used to determine the contribution to signal 
in each sample. The J/V'/o(980) mass distribution shown 
in Fig.rj yields a fitted B° s mass of 5.3748 ± 0.0036 GeV 
and 590 ± 84 B® events, where the uncertainties reflect 
just the statistical uncertainties on the fit. 

Using identical event selections, except for the /o(980) 
mass requirement, a clear J /ip4> peak is found, as shown 
in Fig. [3] The /i + fj," K + K~ mass distribution is fitted 
for a B® — > J/ip4> signal using a double Gaussian func- 
tion with a second-order polynomial for background. An 
unbinned likelihood fit to the J /ip4> distribution shown 
in Fig. E] yields a B° s mass of 5.3631 ± 0.0008 GeV and 
2929 ± 62 B® events, where again the uncertainties are 
statistical only. 

MC signal samples arc used to determine the efficien- 
cies of reconstructing the two B° s decay modes. To take 
account of changes in the instantaneous luminosity, the 
MC samples are overlaid with data events from ran- 
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FIG. 1: The invariant mass distribution of /o(980) candidates 
when the J/ip% + -K~ invariant mass is within ±2 sd of the 
fitted mean B° mass. The solid line represents the fit to all 
the data, and the dashed line the fitted /o(980) signal (see 
text). The vertical dashed lines indicate the region 0.91 GeV 
< M^- < 1.05 GeV. 
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FIG. 2: The invariant mass distribution of B® candidates 
when the tt + tt~ invariant mass is consistent with that of a 
/o(980) meson, i.e., 0.91 < M w+n - < 1.05 GeV. The solid 
line is the fit to all the data and the dashed line the fitted 
signal. The dotted line is a Gaussian function used to 
describe partially reconstructed B decays, (see text). 
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FIG. 3: The invariant J/ipK + K mass distribution when the 
K + K~ invariant mass is consistent with a <f> meson, i.e., 1.01 
< m k+k- < 1-03 GeV. The solid line is the fit to all the data 
and the dashed line is the part fitted to 5° signal, (see text). 
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TABLE I: Relative 


reconstruction efficiencies for different 


running periods. 




Run period 


&reco j Sreco 


Run Ha 


1.19 ± 0.03 


Run Ilbl 


1.29 ± 0.04 


Run IIb2+IIb3 


1.20 ± 0.05 



dom beam crossings collected during each run period. 
In the generation of both the B® — > J/i/j4> an d the 
—> J/^>/o(980) MC signals, a preselection require- 
ment of px > 0.4 GeV is imposed on both kaons and 
pions from the <f> and /o(980). Since the px distributions 
of pions and kaons differ, the preselection efficiencies are 
determined from two additional MC sets of events gen- 
erated without pt cutoffs. 

Reconstruction efficiencies depend on the data-taking 
period (Run Ha - IIb3) as instantaneous luminosity, ag- 
ing of the detector, and changes to the reconstruction 
algorithms affect detector performance. The reconstruc- 
tion efficiencies are therefore measured separately for 
each running period. The instantaneous luminosities for 
data taken during Run IIb3 are similar to those of Run 
IIb2, and the reconstruction efficiencies for Run IIb2 are 
therefore also used for Run IIb3 data. Although the abso- 
lute reconstruction efficiencies depend on the running pe- 
riod, the relative reconstruction efficiencies given in Ta- 
bleware stable. The differences in relative reconstruction 
efficiency are used to estimate a systematic uncertainty 
on The mean relative reconstruction efficiency is 

1.20 ± 0.04, where the uncertainty is from statistics in 
the MC. 

The B° s — > J/tp(f> time development reflects a mix 
of two exponential functions with relative slope values 
driven by the difference in decay widths of the two mass 
eigenstates (Ar,,.). The relative efficiency of any cutoff on 
proper decay length for the two states depends on Ar s 
and the lifetime of the CP-odd eigenstate. The MC sam- 
ples used to determine the relative efficiency use Ar s = 
and the PDG value of T s [T3|. For this Ar s , and assum- 
ing no CP violation, the effect on the relative efficiency of 
/o(98O)/0 is found to be small (~2.5%) and well within 
systematic uncertainties, and therefore no correction is 
applied. 

The branching fraction of B® — > J/ip /o(980) is mea- 
sured relative to B° — > J/ip(p, so any backgrounds that 
peak under the B® — > J/ip(f> mass distribution will affect 
the measurement of R^/^. Possible S- wave contribu- 
tions can arise from the /o or from non-resonant K + J 5 ^ 
production, but these contributions provide only slowly 
varying contributions under the <fi mass peak. The excess 
for larger M K + K - is extrapolated under the <p mas s, g iv- 
ing a possible 5-wave contribution of (12 ± 3)% [24j of 
the total B° a — > J j ip<fi yield. The relative branching ratio 
is therefore scaled up by a factor of 1/0.88 to account for 



an iS-wave contribution to B° s — > J j ijxj). 

One possible background that can affect the observed 
B° s -> J/^>/ (980) yield is the three-body decay B° s 
J / ipTT + tt~ . This background is studied by measuring the 
B° s yield for 7r + 7r~ invariant masses less than the /o(980) 
mass. The 7r + 7r~ mass distribution from non-resonant 
B a s — » J/i/m + ir~ is broad, and measuring the B® yield 
for a sideband in M^+^- therefore provides an extrapola- 
tion of the non-resonant 7r + 7r~ background to the /o(980) 
signal region. In defining a 7r + 7r _ mass window to study 
this background, it is important to avoid regions where 
other known resonances, e.g., B® — » J/ipK*,K* —> Kir 
(with the kaon assigned the pion mass) can contribute. 
The 7r + 7r~ mass window of 0.8-0.9 GeV is chosen be- 
cause this mass range has no overlap with B® — > J/ipK* 
events. 

The mass distribution of /x + /i~7r + 7r~, for 0.8 < M^+ v - 
< 0.9 GeV, is fitted with a floating contribution from 
non-resonant J /ipTT + i:~ decays. The mean and the width 
of the B® peak are constrained to the values obtained 
from the corresponding fit in the /o(980) signal region. 
The fit yields 42 ± 49 events, indicating no evidence of 
Bg — > J/%jm + iT~ non-resonant background, and conse- 
quently no such correction is used in this analysis. 

To check that the results of the analysis do not depend 
on the specific choice of the selection critera, each cut is 
changed around its nominal value, and it is observed that 
Rf /4> does not depend significantly on the exact choice 
of selections. 

The large backgrounds arising from particle combina- 
torics and from partially reconstructed B decays provide 
significant distortion and uncertainties in the distribu- 
tions of background. We study this using same-charge 
pions and the mass distribution from /i + /z~7r ± 7r ± events. 
However, we find that the [i^ pT tx^^ distribution does 
not describe the measured background in our signal sam- 
ple and we therefore do not use it to help constrain the 
distribution of the background. Instead, different param- 
eterizations are used (third-degree polynomial and an ex- 
ponential) to describe the background, and different mass 
regions over which the fit is performed are used to deter- 
mine the signal yield variation. A large variation in the 
number of signal events for B® — > J /ijjf b(980) is found 
for different parameterizations of the background, indi- 
cating that modeling of the background has substantial 
ambiguity. The choice of background parametrization 
comprises the largest contribution to the total system- 
atic uncertainty on Rf a /<t>- 

A similar study of fitting choices is performed on the 
B® —> J/ip4> sample. However, since these backgrounds 
are much smaller and easier to describe, the measured 
event yields change by less than 1%. The presence of 
a B Q — > J/ipir + Tr~ contribution is checked by including 
this channel in the fit, yielding a fit consistent with no 
events. 

The MC distributions of the kinematic variables do not 
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TABLE II: Systematic uncertainties in the branching fraction 
ratio, Rf /4,- 



Source 


Uncertainty 


Fitting 


17.3% 


MC efficiency 


9.2% 


Modeling variables in BDT 


8.9% 


/o(980) mass window 


4.0% 


5-wave contribution 


3.5% 


Total 


22.2% 



describe the data perfectly in all variables. To study this 
effect on the training of the BDT, the MC distributions 
for signal are weighted to match the B® — > J /ipcj) data. 
Only the B a s — > J/ip<j> events are used for this purpose 
because in the B® — > J/tpfo(980) channel there is much 
background and a far smaller signal fraction. 

Using the Run IIb2 data and Run IIb2 MC, we 
find that the relative efficiency for event reconstruction 
changes from 1.20 ± 0.05 without rewcighting to 2.00 
± 0.07 after weighting. Although this corresponds to a 
large difference in relative efficiency, the relative yields 
also change, thereby changing i? /„/</, by just 17.8%. Half 
of the difference between the nominal result and the 
reweighted BDT result is taken as a systematic uncer- 
tainty on Rf /<f,- A 4.0% systematic uncertainty is as- 
signed for the observed dependence of Rf /<p on the size 
of the /o(980) mass window. Table ITTI summarizes the 
values of the systematic uncertainties on R^/^. 

Based on 8 fb _1 of data, DO has extracted a measure- 
ment of the relative branching fraction Rf /,p of Eq. 1. 

R fo/4> = 0.275 ± 0.041 (stat) ± 0.061 (syst). 

This agrees with theoretical expectations and with pre- 
vious measurements of the ratio of widths. 
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